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A new method is proposed for the calculation of the free 
energy of an TV-component <f? 4 theory at finite temperature. 
The method combines a perturbative treatment of the hard 
modes with a non-perturbative treatment in the effectively 
three-dimensional sector of the soft modes. The separation 
between hard and soft modes is achieved by the effective field 
theory method of Braaten and Nieto (ij. One- and two-loop 
gap equations for the screening mass are used to resum higher 
order effects in the soft sector. The proposed method is sim¬ 
ilar to the screened perturbation theory of Karsch, Patkos, 
and Petreczky || albeit avoiding the difficulty of solving gap 
equations in the full four-dimensional theory. Since gap equa¬ 
tions are only used in the three-dimensional effective theory 
the tedious evaluation of Matsubara sums is not necessary. 
This simplification makes it possible to go beyond the one- 
loop gap equation. The large TV as well as the finite TV case 
are discussed. 


In the last few years there has been a lot of progress 
in the calculation of the free energy of high temperature 
field theories @U- On the perturbative level all cor¬ 
rections to the ideal gas up to order g 5 have been calcu¬ 
lated for gauge theories 11 as well as for scalar <T> 4 the¬ 
ory PP1 . Unfortunately, these calculations revealed 
that already for moderate values of the coupling con¬ 
stant higher order corrections become large and oscillat¬ 
ing, i.e. each additional higher order correction flips the 
sign of the deviation from the ideal gas limit (cf. fig. j^). 
In addition, the results depend crucially on the chosen 
renormalization point j^|| . 

In this talk I would like to contribute to the question 
whether/how the perturbative calculation of the free en¬ 
ergy can be improved such that one gets a reliable an¬ 
swer not only for very small values of the coupling con¬ 
stant (where perturbation theory works, but there one 
might use the free (unperturbed) result as well), but also 
for moderate values of g. The aim is to reorganize the 
loop expansion such that higher loop corrections do not 
drastically change the result. Of course, it is not clear 
whether such a reorganization is possible or if one has to 
rely solely on genuine non-perturbative approaches like 
lattice calculations (24). Even for theories where lattice 
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calculations yield reliable results an improved perturba¬ 
tive calculation scheme is nonetheless useful to identify 
the relevant degrees of freedom of the system. Indeed, 
the analysis of QCD lattice calculations indicates that at 
high temperatures the system of massless coupled modes 
can be very well approximated by a gas of noninteracting 
massive particles |^5| ^8j. It is clearly desirable to sup¬ 
port this heuristic finding by a calculation starting from 
first principles such that also the fit parameters (e.g. the 
temperature dependence of the mass) can be calculated 
from the underlying theory. In this spirit a reorganiza¬ 
tion of the loop expansion called screened perturbation 
theory was suggested in ( 2 ) for an TV-component scalar 
4> 4 theory. 
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FIG. 1. The perturbative result for the free energy density 
of a one-component massless g 2 ® 4 theory as a function of the 
coupling g. The label n marks the highest perturbative con¬ 
tribution ~ g n which is taken into account. The free energy 
density is normalized to its ideal gas value. 

Starting from the Lagrangian 

£ = i ( 5^ )2 + ^ ( ^ )2 ( 1 ) 

(where i labels the TV components of the fields) reorgani¬ 
zation of perturbation theory is achieved by adding and 
subtracting a mass term: 



G — Un + C, 


int 1 


C 0 = l -{d^) 2 + l -M^l 
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The mass M is calculated from the one-loop gap equation 


M = 


Q 


(3) 


In || the necessary loop calculations were performed up 
to three loops for the large N case and up to two loops 
for arbitrary N. The result for the free energy density 
for the large N case is shown in fig. || taken from 0. 
Obviously the convergence of the screened perturbation 
theory looks very promising even for large values of the 
coupling constant. 



FIG. 2. The free energy density of the scalar held theory 
in the large N limit as a function of the coupling constant in 
units of the free energy density of a massless ideal gas. The 
curves represent three-loop (a), two-loop (b) and one-loop (c) 
results of the screened loop expansion as well as the 0(g i ) 
result (d) of the conventional perturbative expansion. For 
details see [|2j] where the figure is taken from. 

Of course, it is interesting to check whether this scheme 
works equally well for arbitrary N. There is a techni¬ 
cal reason why the three-loop calculations were not per¬ 
formed in H for small N: It is very complicated to cal¬ 
culate the sunset and the basketball diagram shown in 
fig- 1 a and b, respectively. While the latter contributes 
to the free energy, the former shows up as a two-loop 
correction to the gap equation !)• In the large N limit 
both diagrams are suppressed by a factor l/N. For finite 
N, however, these diagrams have to be taken into ac¬ 
count. Unfortunately, to calculate these diagrams with 
massive propagators is quite messy especially to perform 
the Matsubara sums. 

To circumvent this technical problem the following ob¬ 
servation from the QCD case might be useful: If the (os¬ 
cillating) perturbative corrections to the massless ideal 
gas are separated in their contributions from the hard 


(order of temperature T) and from the soft (order gT) 
sector H one observes that the mentioned oscillatory be¬ 
havior is mainly caused by the soft modes, i.e. for not too 
large values of the coupling constant the contribution of 
the hard modes to the free energy is well behaved while 
the one from the soft modes is not. The latter, however, 
can be described by an effective three-dimensional (vac¬ 
uum like) theory. Hence, if non-perturbative improve¬ 
ments are restricted to the soft mode sector, Matsubara 
sums do not show up and one can use the techniques of 
vacuum perturbation theory to calculate the necessary 
integrals which enter the respective gap equation and the 
free energy calculation. This consideration is the basis of 
the work presented here. 
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FIG. 3. Sunset (a) and basketball (b) diagram contribut¬ 
ing to the gap equation for the mass and to the free energy 
density, respectively. 

The recipe to calculate the free energy density is the 
following: 

1. Disentangle the contributions from different scales 
(hard and soft). 

2. Use conventional perturbation theory to calculate the 
contribution from the hard scale. 

3. Use screened perturbation theory to calculate the con¬ 
tribution from the soft scale. 

Fortunately, the first two points have already been 
performed for one-conrponent <f> 4 theory in 0 and for 
gauge theories in p9| . For simplicity and to compare my 
results with the four-dimensional screened perturbation 
theory developed in 0 I will work in the following with 
an iV-component <f> 4 theory which is a straightforward 
generalization of the system studied in |jl[]. The free en¬ 
ergy density receives contributions from the hard modes, 
Fh , as well as from the soft modes, F s . In conventional 
perturbation theory up to 0(g 5 ) these contributions are 
schematically given by[J 

F h = # + #g 2 + #<? 4 (4) 

F s = #ff 3 + #ff 4 + #S 5 (5) 

The perturbative contribution from the hard sector is 
given by (cf. 0) 

F h = F?+Ff + F% 1 ( 6 ) 



with 


1 Strictly speaking the number sign in front of g 5 represents 
not a pure number but contains also a log g contribution. 
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where g denotes the renormalization scale, N= (7V+2)/3, 
7 is Euler’s constant, and £( 2 ) is the Riemann zeta func¬ 
tion. 

The three terms contributing to F s in (|s]) come from 
one-, two-, and three-loop diagrams calculated in the ef¬ 
fective three-dimensional theory [Q 
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The so-called short distance coefficients m and A in this 
theory for the soft modes are influenced by the interac¬ 
tion with the hard modes. Perturbatively they are given 

by 


A = g 2 T + o(g 4 ) 


( 11 ) 
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where the higher order corrections do not influence the 
result for the free energy up to 0(g 5 ). Here a new pa¬ 
rameter shows up, the separation scale A which separates 
the hard from the soft modes. A serves also as the renor¬ 
malization scale in the three-dimensional theory, i.e. it 
appears in the course of renormalizing logarithmic ul¬ 
traviolet divergences from loop integrals. Note that the 
short distance coefficient m is not a physical quantity. If 
physical quantities (like the free energy or the screening 
mass) are calculated in conventional perturbation theory 
the dependence of the short distance coefficients on A ex¬ 
actly cancels in every order of g the A dependence com¬ 
ing from the renormalization of loops |jj . Going beyond 
a perturbative calculation in the following a dependence 
on A will remain. This will be discussed below. 

According to point 2 of my recipe I will keep the pertur¬ 
bative results caused by the hard modes, i.e. the results 
for Fh, m and A, while according to point 3 I will replace 


the three terms of F s given in (||) by the one-, two-, and 
three-loop result of (three-dimensional) screened pertur¬ 
bation theory 


£ — A + Ant ) 

A = \(d<t>i ) 2 + \M 2 <t> 2 i , (13) 

Ant = -\( M2 - ™ 2 )^i + ^V^i ) 2 • 

The Feynman diagrams which have to be calculated are 
depicted in fig. [|. The cross denotes the counter term 

5m 2 = M 2 - m 2 (14) 

(times —1). There is an additional counter diagram not 
shown in fig. which serves to renormalize the mass. 
It cancels the singularity arising from the basketball dia¬ 
gram. It also yields a finite contribution to the three-loop 
result given below. For details of the calculation cf. [ j30| . 
The result for the contribution of the (resummed) soft 
modes to the free energy density is 

F s = T{fl l + f 2l + ff) (15) 

with 

fl l = -^-NM\ (16) 
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where the appearance of a cross in a diagram is treated as 
being a loop, i.e. the third diagram of fig. ^ is a two-loop 
contribution and the fifth and sixth diagram are three- 
loop contributions. 
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FIG. 4. Feynman diagrams for the free energy contribu¬ 
tions to be calculated in three-dimensional screened pertur¬ 
bation theory with massive propagators with mass M. The 
cross denotes a counter term —5m 2 . 


The final result for the free energy density is 
F = F h + F„ 


(19) 


where Fh is given by (pN9|) and F s by (15-lq) with the 


parameters given in ( |T7| . |I 2| . |1 l| i . By inspecting these equa¬ 
tions we find that F depends on the temperature T, the 
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coupling constant g , the renormalization scale (of the full 
four-dimensional theory) g, the separation scale A, and 
on the so far unspecified mass M. In a full calculation of 
the free energy density only the dependence on T and g 
would remain. In an approximate calculation, however, 
there is an additional dependence on the other three pa¬ 
rameters and appropriate values have to be chosen for 
them. Since I am mainly interested in the soft mode sec¬ 
tor and in a proper choice for M I will determine g and A 
in a very simple way: I require that perturbation theory 
in the hard mode sector works optimal, i.e. I will choose g 
and A such that the g 4 contributions to F/,, given in (|fj), 
and 77i, given in (12), vanish. This fixes the separation 
scale to 


A = 4»re Xp| -“_AAi! + I<V!l 
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while the renormalization scale becomes 
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This yields e.g. 


zeroth and the first/second order. Since the free inverse 
propagator vanishes at this momentum this simply means 
that the self energy is demanded to vanish at k 2 = — M 2 . 
Taking the second order result for the self energy, i.e. up 
to two loops, corresponds to the evaluation of the free 
energy up to three loops. This yields the gap equation 
(see also |il] [h| and references therein) 
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fc 2 = -M 2 • 

(25) 


If the calculation of the free energy is restricted to max¬ 
imal two loops the corresponding one-loop gap equation 
results from ( f25[) by simply dropping the last three di¬ 
agrams. Again, there is an additional counter term not 
shown in ( |l5| ) which renormalizes the mass. It cancels the 
singularity arising from the sunset diagram. For details 
of the calculation see Q. The result for the two-loop 
gap equation (TLG) is 
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(TLG). (26) 


gpz 1.12 • 27tT for TV —> oo (22) 

and 

g « 0.57 T for TV = 1. (23) 

Admittedly, the chosen value for g is somewhat low for 
the case TV = 1 since it is usually supposed to be some¬ 
where around 2nT (cf. HII-FRl )■ However, changing g 
or A does not qualitatively change the results to be pre¬ 
sented in the following. 

Of course, the crucial question now is how to deter¬ 
mine M such that large perturbative contributions from 
higher loops are resummed in the tree level mass M. I 
will explore two different resummation procedures: 

The first one is based on the principle of minimal sen¬ 
sitivity (see e.g. and references therein). Since the 
full calculation for the free energy is independent of M 
it might be reasonable to look for an extremum of the 
approximate result with respect to M, i.e. 

A = 0 (POMS). (24) 

The second resummation procedure is based on the crite¬ 
rion of fastest apparent convergence (see e.g. |3l]] and ref¬ 
erences therein). Here, an appropriate quantity is chosen, 
calculated in two different orders of the approximation 
scheme, and it is required that the difference vanishes. I 
will use the inverse propagator evaluated at momentum 
k 2 = —M 2 and take the difference of the results for the 


For the one-loop gap equation (OLG) all terms on the 
right hand side except the first two have to be dropped. 

Both schemes to determine M resulting in (^4j) and 
([2(f), respectively, might yield more than one solution for 
M. I always take the one which is continuously connected 
to the perturbative result ( M 2 « to 2 ). If the calculation 
of the free energy is restricted to maximal two loops the 
principle of minimal sensitivity yields the same result as 
the one-loop gap equation. In the large TV limit the sunset 
diagram (the last one in (|25|)) is suppressed. There TLG, 
OLG, and POMS yield identical results. This is easy to 
understand, if one recalls that for TV —> oo the tadpole 
resummation (which is already achieved by OLG) is suf¬ 
ficient to solve the exact Dyson-Schwinger equation for 
the propagator |l7). If the solution for M is inserted in 
SI) the three-loop contribution ff exactly vanishes 
in the large TV limit. The same is true for all higher loop 
contributions not calculated here (see also 0 ®)- 

For the presentation of the results one has to make a 
decision how to sort the various contributions Ff, F 21 , 
Ff ■ fs l i ff i and fs l • O ne might be tempted to sort them 
in numbers of loops such that the lowest order approxi¬ 
mation to the free energy density would be F t \ l + T/ S i; , 
the next one F^ 1 + Tf] 1 + F 21 +Tf 21 , and so on. This, 
however, I think is misleading in view of the schematic 
picture of the perturbative contributions given in (|JJ^) 
and the fact that the contributions / s 1( , f 21 , and ff only 
improve the soft perturbative g 3 , g 4 , and g 5 contribu¬ 
tions, respectively. Therefore, I choose the following se¬ 
quence of approximations: 
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(31) 


Note that for the results shown below the renormalization 
scale /i is chosen such that F^ 1 vanishes. Therefore the 
difference between F 3 and F 4 is solely due to ff. 

A second question which has to be clarified is which M 
is taken for which order of the approximation. Concern¬ 
ing the principle of minimal sensitivity for the approx¬ 
imation F 3 , clearly, the mass M has to be determined 
from (|24l) with the input F = F$. With the resulting 
mass, F 5 can be evaluated as a function of the coupling 
constant g (and the temperature T). One might take the 
same mass for the evaluation of F4. However, I think it 
is more appropriate here to use only F = F4 as an in¬ 
put in (|24). This philosophy, however, does not work for 
F 3 . Since fj l is proportional to M 3 , POMS would yield 
M = 0, i.e. no difference between the approximations F 2 
and F 3 . This is, of course, not very useful. Therefore I 
break the rule for this special case and use for the eval¬ 
uation of F 3 the value for M resulting from (|24|) with 
F = F4 as input. 
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FIG. 5. Various approximations F n (labeled with n) to the 
free energy density of an A-component 4 > 4 theory as a func¬ 
tion of the coupling g in the large N limit. See main text for 
details. The free energy density is normalized to its ideal gas 
value Fq . The dashed lines show the perturbative approxima¬ 
tions 0(g 3 ), 0(g 5 ), 0(g 4 ) (top to bottom). 

For the criterion of fastest apparent convergence and 
its resulting one- and two-loop gap equations I will use 
the result of TLG for the evaluation of F$ and of OLG 
for the evaluation of F4. Again, it is not clear which M 
should be used for F 3 . I will use the result of OLG here. 

The results for F 0 - F 5 are shown in fig. g for the case 
N —> 00 and in figs. g and g for TV = 1. To make the 


two cases comparable the respective maximal value of 
g is chosen such that the short distance coefficient m 
given in © is approximately the same for both cases. 
For the large TV case I can reproduce the results of || 
(cf. fig. g): Perturbation theory is improved by the re¬ 
summation scheme; even for large values of the coupling 
constant the best approximation (full line in fig. ! f 5 
in fig. g) is quite near to the ideal gas limit (about 10 % 
deviation for g = 12); the screening mass M (not shown) 
equals the perturbative one-loop result for small coupling 
constants (g 2 part of the short distance coefficient (|l^)) 
and is about half of it for g = 12 . 
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FIG. 6 . Various approximations F n (labeled with n) to the 
free energy density of a one-component l ? 4 theory as a func¬ 
tion of the coupling g using the criterion of fastest apparent 
convergence. See main text for details. The free energy den¬ 
sity is normalized to its ideal gas value Fo. “pert” labels the 
perturbative approximation 0(g 5 ). 

Unfortunately, the picture is quite different for the 
case TV = 1. Comparing the purely perturbative re¬ 
sults (fig. g) with the ones obtained in the resumma¬ 
tion schemes (figs. !@) it turns out that the criterion of 
fastest apparent convergence (fig. g) does not improve the 
loop expansion at all. The principle of minimal sensitiv¬ 
ity improves the loop expansion. However, it is striking 
that also in this scheme for larger values of the coupling 
constant the approximate value for the free energy den¬ 
sity changes from \F\ < |F 0 | to |F| > |F 0 | when pro¬ 
ceeding from the fourth to the fifth approximation. This 
seems to resemble the oscillatory behavior of the per¬ 
turbative results. Obviously, there is a window between 
g ss 3 and g ss 5 where the latter resummation scheme 
(POMS) shows a promising convergence behavior and im¬ 
proves the perturbative results. (Below g ss 3 one might 
trust the perturbative results as well.) However, the fact 
that the two discussed resummation schemes yield quite 
different results raises doubts on the reliability of both. 
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FIG. 7. Same as fig. |g| but using the principle of minimal 
sensitivity. 

Before discussing the possible lessons one could learn 
from the three-dimensional version of resummed pertur¬ 
bation theory let me recall the additional assumptions 
which enter the presented method besides the general 
idea of resumming the soft modes only: First of all, I had 
to choose the values for the renormalization scale g and 
the separation scale A. I have used a very simple recipe 
here to fix these quantities by demanding that the g A con¬ 
tributions to the hard mode quantities F h and m have to 
vanish (cf. (p(i|-^3|) and discussion there). It is worth not¬ 
ing that this approach is in accordance with the criterion 
of fastest apparent convergence which I have used later 
on. In principle, if perturbation theory works these g 4 
contributions would be small but in general would not 
completely vanish. Therefore, forcing them to vanish by 
choosing the scales accordingly might cause other contri¬ 
butions to become unnaturally large. On the other hand, 
this simple choice allows to discuss in the cleanest way 
the convergence behavior of the soft mode sector and its 
influence on the total result for the free energy. If e.g. g 
is chosen in another way then the influence of f 21 and 
ff 1 in (|3(^,|3l]) would intertwine with the influence of F% 1 . 
Nonetheless, it is important to check whether at least the 
qualitative results are robust against changes in g and A. 
Therefore, I will now explore other choices for the scale 
parameters. I will restrict myself to the most interesting 
case N = 1. 

Concerning the renormalization scale a reasonable 
choice might be g = 2irT since this is the typical en¬ 
ergy scale of the hard modes @I0- For this case the 
contribution of the hard modes to the free energy (b[) 
becomes (for N = 1) 

F h « Ff [1 - 0.0079 g 2 (l - 0.046 g 2 )] for g = 2irT. 

(32) 

One finds that the g 4 contribution overwhelms the g 2 
contribution already for g ss 4.7. In view of the large 


coupling constants I have dealt with in this work (up to 
g = 7 for N = 1) this choice of g is inappropriate since in 
this case the naive perturbation theory which I have used 
in the hard mode sector breaks down. For the more naive 
choice g = T one gets a g 4 coefficient which is somewhat 
smaller: 

F h « F^ 1 [l — 0.0079 g 2 (l - 0.011 g 2 )] for g = T. 

(33) 

For this choice of g (and A as given in (^0|)) the free en¬ 
ergy contributions are plotted in fig. using the principle 
of minimal sensitivity which for the original choice of g 
was the better one as compared to the criterion of ap¬ 
parent convergence. Obviously, also for this choice of the 
renormalization scale POMS improves the convergence 
as compared to the perturbative loop expansion. How¬ 
ever, one finds qualitatively the same feature as in fig. [?], 
namely the oscillatory behavior around the ideal gas limit 
for coupling constants larger than g ss 5. Note, however, 
that the fourth and fifth approximation basically have 
now exchanged their places as compared to fig. [t|. 
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FIG. 8. Same as fig. ^ but using the principle of minimal 
sensitivity and g = T. 

Concerning the separation scale it is reasonable to put 
it somewhere between the typical energy scales of the 
hard and the soft modes. While the former is 2nT as dis¬ 
cussed above the latter is characterized by m or roughly 
the lowest perturbative ( g 2 ) contribution to m (as long 
as perturbation theory holds for the hard mode sector): 

m U : = ^l^ 92 (F)T 2 ■ (34) 

Thus a reasonable choice is 

A = \J2-kT mu ~ y/g . (35) 

Fig. H shows the free energy contributions for this choice 
of the separation scale (and g as given in (p3|)) using 
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again the principle of minimal sensitivity. Comparing 
the values for F 5 from fig. || and fig p| shows that the 
latter choice for the separation scale ( |3fi| ) makes the con¬ 
vergence behavior worse. In addition, it turns out that 
in this case the full POMS condition, i.e. using F = F$ 
in (HJ) , does not have a solution for M for coupling con¬ 
stants smaller than g ss 2. Therefore, in fig. F 5 is plot¬ 
ted only for g > 2 (hardly visible). All of this is qualita¬ 
tively in line with the general finding that the presented 
resummation scheme does not improve the loop expan¬ 
sion for the case N = 1. 
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FIG. 9. Same as fi g. |(i| but using the principle of minimal 
sensitivity and A = \Z2ivT mu- 

Besides the specific choices for the renormalization and 
separation scales I have tacitly assumed that there is only 
one counter term ( |l4| ) in (|l^). Alternatively one could 
introduce an additional counter term for each additional 
loo p or der. This approach will be discussed in more detail 
in Q. Qualitatively this modified method yields the 
same results as presented here. Especially no further 
improvement concerning the poor convergence behavior 
for N = 1 is achieved by this modification. 

To summarize, I have found that a scheme which 
combines conventional perturbation theory for the hard 
modes with screened perturbation theory for the soft 
modes significantly improves conventional perturbation 
theory for an IV-component 3> 4 theory in the large N 
limit. Contrary, this scheme does not work equally well 
for finite N. There are several possibilities to explain this 
shortcoming: 

1. The presented scheme is based on the assump¬ 
tion that the hard and the soft modes can be sepa¬ 
rated and that the screening effect can be neglected for 
the hard modes. The necessary condition for this as¬ 
sumption to hold is that in the Matsubara propagator 
1 /{k 2 + (2 mrT) 2 + M 2 ) the mass term can be neglected 
compared to the temperature term for n /0, i.e. 

M 2 < (2ttT) 2 . (36) 


If the screening mass does not fulfill this condition 
the presented scheme would not be appropriate. In 
turn, if screened perturbation theory works in the four- 
dimensional theory and yields a screening mass which 
fulfills (|3^) then also my three-dimensional version of 
screened perturbation theory should work. Of course, 
I cannot say what the true value for the screening mass 
is. However, I can at least check whether the result for M 
obtained in my approximation scheme is consistent with 
the underlying assumption (|3(]). Fig. ^o] shows the re¬ 
sults for M for the different resummation schemes. The 
line labeled with “POMS” is obtained from ( |24| ) using 
the “best” result F = F 5 from (fH|). The “TLG” and 
“OLG” lines result from (|2^) with the appropriate terms 
dropped in the latter case. I think it is fair to say that the 
POMS and OLG results fulfill ( ^6| ) in the plotted regime 
for the coupling constant. For the TLG result relation 
( p6| ) is not very well satisfied for the largest plotted values 
of the coupling constant. However, for g « 5 there is a 
factor of about 40 between M 2 and (27 tT) 2 so that ( |36| ) is 
roughly fulfilled. Fig. || shows that already here the con¬ 
vergence of the resummation scheme is worse than the 
conventional perturbative scheme. Thus, the failure of 
the resummation scheme utilizing the criterion of fastest 
apparent convergence cannot be traced back to its incon¬ 
sistency with the underlying assumption 

2. Another reason for the finding that the presented 
resummation scheme does not work equally well for N —> 
00 and N = 1 may be that only a momentum indepen¬ 
dent mass is resummed. It might appear that the re¬ 
summation of a more general, i.e. momentum dependent 
self energy (as e.g. advocated in Q) would be more ap¬ 
propriate. For the large N case this would not change 
anything since the sunrise diagram and all higher mo¬ 
mentum dependent diagrams are suppressed by powers 
of 1/N. As already pointed out before, the resummation 
of tadpoles is the only thing which has to be done in the 
large N limit. If the possible momentum dependence is 
a necessary ingredient for a reliable resummation scheme 
then this would naturally explain why the scheme pre¬ 
sented here works for N —» 00 but fails for N = 1. I will 
elaborate on this question in Q. 

3. Besides its momentum dependence the sunrise di¬ 
agram (and the corresponding basketball diagram con¬ 
tributing to the free energy) gives rise to a logarithmic 
term in the mass M. One may argue that this loga¬ 
rithmic term via an inappropriately large contribution 
ln(#A/M) shifts the mass to a wrong place and causes 
a breakdown of the three-dimensional screened pertur¬ 
bation theory. Since this logarithmic term is caused by 
the renormalization of the theory it might be possible to 
handle that problem by renormalization group equations 
(cf. e.g. §). Of course, this is a pure speculation so far 
and deserves further studies. If it is true it would also 
explain why the resummation works in the large N limit. 
There, these logarithmic contributions are suppressed as 
one can easily check by inspecting (|Ii|) and (|26|). 
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FIG. 10. The resummed mass M as obtained from POMS, 
TLG, and OLG and the short distance coefficient m as a 
function of the coupling g for a one-component <E> 4 theory. 
See main text for details. 

4. The most negative explanation for the shortcomings 
of the presented resummation scheme would be that the 
concept of screened perturbation theory does not work at 
all (or only for very special cases like the large N limit) or 
more generally that concerning the determination of the 
free energy there is so far no method which reorganizes 
the loop expansion in a useful way, at least for coupling 
constants which are not tiny. This would imply that we 
have to rely on genuine non-perturbative approaches like 
lattice calculations. Indeed, the free energy density is 
a quantity which can be calculated on the lattice (see 
e.g. |||,^|). However, for many other, e.g. non-static 
quantities there are fundamental obstacles which so far 
prevent their evaluation on the lattice. If (improved) per¬ 
turbation theory fails to calculate the free energy density 
one may also doubt the results for other quantities ob¬ 
tained in a loop expansion in high temperature quantum 
field theory. 

I conclude that the question how to reliably calculate 
the free energy density in a modified loop expansion is 
still open. 
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